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(iii) 


| NTRODUCT | ON 


The modern theory of continuum mechanics has expanded greatly in 
the last few years. As interest in finite deformation theory increased, 
mathematical fundamentals for general types of materials were established, 
and the techniques of differential geometry were brought to bear on the 


problems. 


In this thesis, we outline some of the general theories, and give 
a specific application. We generalize a problem of Vaughan's [5], and dis- 
cuss the deformation of a compressed, laminated cylinder under its own 


weight. 


In the first part, we rely on the concepts introduced by Wang [6], 
[7] to give a mathematical description of laminated bodies. We develop the 
theory in terms of local coordinates, and confine our discussion to simple, 


elastic, materially uniform, isotropic bodies. 


In the second part, we use the techniques of Green and Zerna [2] 
to deal with the problem of smal! deformations on large. We develop an 
appropriate notation for the case of laminated or pre-stressed bodies. After 
deriving the general equilibrium equations for a perturbation, we snow how to 
apply the theory to the stress-system of a universal solution. We choose a 
particular constitutive equation and a particular metric, to simplify The 


resulting equations. 


Finally, we use numerical techniques to obtain solutions to the 
equations, subject to appropriate boundary conditions. We show that lamina- 


tions have a considerable effect on the deformation. 
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(iv) 


The summation convention is used throughout this thesis, unless 
specified otherwise. Although we have tried to use a consistent notation, 


the familiar notation for classical quantities has been retained. 
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LIST OF SYMBOLS 


a,b,u,v - elements of the isotropy group. 

A, - the values of the coordinates x, Ofraspaini= a in G . 
B - the left-Cauchy-Green tensor. 

By - the left-Cauchy-Green tensor of the lamination. 

ay FBF! 


B - an adjusted form of B , defined by (2.2.35). 


B. - one of the vectors spanning Hs 

B - the physical body. 

C - the right-Cauchy-Green tensor. 

e578, u - linear frames at p in M., 

ey 6P) - the linear frame corresponding to the standard frame of 
EB (a) dbase Vector ar a in G . 

Fy (2) - a base vector at @ in E(M) 

E(M) or E(B) - the bundle of linear frames on M. 

E(M,U) or E(U) - the bundle of reference frames. 


E~(B) or To (E(B) ) -~ the tangent space of E(B) at u. 
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Hs - the horizontal subspace of Ex(B) 


I - the identity tensor. 
I,,1,,1, - the three principal invariants of B. 
I, - in Chapter 3, the modified Bessel function. 
L - the cube of the compression ratio. 
Le - one of the laminae. 
M - a lamination parameter defined by (3.6.1) and (3.6.3). 
M - the lamination parameter for the first harmonic. 
M - the body manifold. 
oF or Nate - the tangent space of M at p. 
p,q. —- points fn °M* 
r_,r_(p) - local configurations of M 
p a Pp 
= - the material isomorphism between M_ and M 
Pq Pp g 
3 : 
R - Euclidean 3-space. 
s - aradial variable defined by (3.4.3). 
ee the boundary value of s_ for the first harmonic. 
, 5 
SL(3) - the special |inear group for R . 


S0(3) - the special orthogonal group for R° : 
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(viii) 


t - the stress vector. 


T - a function of the radial variable, defined by (3.2.1). 
T(M) - the tangent bundle on M. 
T(E(B)) - the tangent bundle on E(B) 
cat ali a base vector for oe at u. 
(gy Fay) - a reference chart on M 
Gy, 11.) & era material chart on T(M) 
U - areference atlas for M . 
V += ascalar multiple of the function T . 
ait - the vertical subspace of E~(B) : 
we u,v - the components of the perturbed disp! acement vector. 
W - the strain-energy function. 
x,y -. vector fields on E(M) 
Yereweacvector, field on G . 
X,Y XX, - Cartesian coordinate systems. 
7 - a function of the vertical variable, defined by C3 a2 oe 
€ - a perturbation parameter. 


« - a lamination parameter defined by (3.6.1). 
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(ix) 


No Tap - bundle maps for T(M) 

n(u) - the material atlas. 

{60},{6% - curvilinear coordinate systems. 

{e~} - convected coordinate system. 

8,6' - atlases for T(M) 

138 - components of the "perturbed-stress" tensor, defined by (2.2.15). 
Ey ~ @ bundle map for E(M) . 

Ops—) density. 

Oo - f divided by the initial height of the cylinder. 

ore - the boundary value of o. 

T - the Cauchy stress tensor. 

$,W - diffeomorphisms of M to R°. 

®,¥ - Mooney-Rivlin coefficients. 

Special Symbols 

Z\5 ~- the natural lift of aie to the frame u. 

“Ap. - the special lift of zj, as defined by (1.5.49). 

(Ly) |, or Lz,y J], - the Lie derivative of y with respect to z. 
VEY - the covariant derivative of y with respect to z. 


Q=Q+t <0 - the perturbed general quantity Q. 
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The Material Geometry of Stress 


1.1 Materially Uniform Simple Bodies. 


The points of a physical body B will! be regarded as points of 
a closed, connected, 3-dimensional, differentiable, orientable manifold 
M . M_ can be covered by one coordinate neighbcurhood, i.e. there exist 


diffeomorphisms 


¢:M>R Cletet) 


where R is Euclidean 3~space with standard orientation. 


The tangent space at p«M_ is denoted by us or TM) ; 


where T(M) is the usual tangent bundle. A map 


is called a local configuration if it is an orientation preserving isomor- 


phism. Each diffeomorphism (1.1.1) induces a map at p 


x 
SAMOS he : es ee 
x|> p 
@ is called a configuration if |, is.a local configuration for all p . 
There need not be any (global) configuration p:M-> R° such that 
Sols. tOree Dahir Lula raat, 
Py p pD p p 


We refer to the manifold as M or B if no confusion will arise. 
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A motion of B is amap @: RxM-> Rt defined by 
O(s,M) = (so. (M)) Clot. @) 


where > (M) is a configuration ¥s . We represent Oe by its deformation 
gradient 
Ee Rox aie 5: F_ = ° ; 
j Ps |p ° Tp 
where - is some local reference configuration. Each choice of - gives 


a different representation of F ... If o! are the coordinates under 6 


: i 
and Bon are the base vectors under rs , then (CFD) = eal A ps a8, 


ar or 
p p 


the chain rule. Since " and %. are orientation preserving, detF>0O. 


We have the following diagram. 
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We define a second-order symmetric stress tensor Tt by 


Oe tls Cot. &) 


where n is the normal to a surface at a poln? -p of Bi, and SS is the 


stress vector at that point. 


B is called a simple body if t is determined solely by the 


local history: 


where t,s denote time and L is the response functional. For any choice 
of reference configuration - at p , we can refer to the deformation 


gradient and write 


t 
T(t pde= oak (Fo,p) : CRESS) 


=—cO 


If the body is elastic, there is no dependence on the history of 


the motion, so we may write 


T(t,p) = T(p) = RF) Gites 


We now develop a geometry for B which is characteristic of its 
response to stress. We say that two points p,q are materially isomorphic 


if there exist local configurations rhe such that 


Rott = ROCF) fot al Fs Ci 10) 
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We take B to be a simple, materially uniform body, where by 


materially uniform we mean that any two points p,q of M = are material 


isomorphic. 


A set U = LOG) ; a e A} is a reference atlas for M if it 


has the following properties: 


{uy : ae At is an open cover of M 


tne > @« A} is a smooth field of reference configurations, 


Py CPi: M > R {OG epeen U 


it ,p;decsuy .G then r7!(q) o r (p) is a material isomorphism. 
fe} a fe} 


bt oe ea ee Ug » then rp (q) © ry fP? is a material isomorphism. 
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Us is maximal. 
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Qn the open set uy we may write 
T(uy) = R, K Fe) ; Gl. be to) 


since by property (3) there is no dependence on the individual points. In 


fact, by property (4), the charts are compatible on the overlaps, i.e. 
UES ea tG,) 
so we may write 


T = RF) : Pt 4) 


We deal exclusiveiy with smooth, simple, materially uniform bodies, 


which may be equipped with a reference atlas. 


1.2 The Isotropy Group. 


The isotropy group 35 is the group of all material isomorphisms 


of p with itself, i.e. 
G=Hh- 6b BB eC a he o h)} é Cl219 
t p 5 | p'ox| px] 5 


In addition, we ask that Seas Shee? =, det h =\+1 


Rather than deal directly with Bo , we can transfer this concept 
to R . We say G: R ee R° is a material isomorphism relative to ro 
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Fig. 1.2.3. 


he 35 . Clearly the stress tensor is intrinsic. 


A tensor field wu is defined to be intrinsic if the element u 
on Se is intrinsic in the above sense and for any two points p,q , we 


have 
Liat Co eect Ti - Clee. }O) 


where SM is a mapping on the tensor space, induced by the material 
isomorphism Mog . Clearly, if H, is intrinsic, we can create an intrin- 
sic field wu through the mappings Pye . Since these mappings are assumed 


Sroary, SO lS — iL. 


lf Cu_,r.) and (u,,r,) are two charts of U on a neighbour- 
a? oO Boas 
hood of p*, then ry SP) ° rg) is a material isomorphism. Define 
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The C48 form the group of coordinate transformations. If =. are the 
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A point p is a solid point if there exists a local configura- 
tion rh > _ ¢ SO(3) , otherwise p is a fluid point. In a materially 
uniform body B , all points are solid if one is, so B is SOU att 
Gee Ue). 1S PSOTLODLG It -o0\5) G.G". Hence for a solid Isotropic 


body, G = S0(3) 


We make one further assumption, necessary for the development 
of the material geometry: the isotropy group is a closed Lie subgroup of 
SL(B,) . This assumption is satisfied by isotropic solid bodies in parti- 


cular. 


1.3 The Specialized Constitutive Equation. 


We now show how the constitutive equation (|1.1.17) can be 
reduced to a simple coordinate form. By the polar decomposition theorem, 


a matrix F , det F > 0, may be factored uniquely: 
Pease => Vege Gleoet) 


where R is orthogonal and U,V are positive definite and symmetric. We 
define two tensors B and C,, called the left and right Cauchy-Green 


tensors, respectively, by: 
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Now for a solid isotropic body, G = SO(3) , so if we choose 


G= R! Intl. 4.9). We gat 


T SURUER PEER CVRR = ROR H(B) |. (1.323) 


In preliminaries to their major work on isotropic functions, 
Rivlin [3] and Ericksen and Rivlin LI] used physical principles (cf. Appen- 


dix A) to show that the actual form of (1.3.3) for an isotropic elastic solid 


must be 


. Z 
T = fot + FB + FB " C5324) 


where the f's are functions of the three principal invariants of B. 


We also state now, and discuss more fully in Section |.6, that 
on a solid body we have an induced, intrinsic, Riemannian metric g which, 


relative to some local reference configuration i » has the form 
ees oe Oe, (1.3.5) 


Here the fact that r may possibly not be induced from a global configura- 


tion is important. 


We assume that we have a global configuration $ which gives an 
X - system of coordinates. We denote the deformation configuration by yp, 
giving an x - system of coordinates, and we also have base vectors a 
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and no reference to r_ appears. 


Now that we have an explicit representation for the stress tensor, 
we need the equilibrium equations. To derive these, we have to discuss 


covariant differentiation on a manifold. 


1.4 Material Bundle Spaces. 


Material bundle spaces are new differentiable manifolds over M 
which are determined by the response of the body B. We will describe their 


construction by reviewing the construction of the ordinary bundle spaces. 


To discuss differentiation on M , we require two new manifolds: 


the tangent bundle, T(M) , and the bundle of linear frames, -G(M) . 


Briefly, a material tangent bundle is a sub-bundle of T(M) , 
characterized by the fact that its structure group is not GL(3) but G, 


the isotropy group. 


The restriction is often significant. Consider a spherical shel] 
of a transversely isotropic material whose axis of transverse isotropy is 
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is a material isomorphism for all p« Uy» IE Ug - A material atlas is 
a maximal collection of pairwise-compatible material charts 


{Cuyny) : a ¢€ A} , where the set {uy : a € A} is an open cover of M. 
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the material atlas. If G is the isotropy group corresponding to U, 


then G is the structure group of the material tangent bundle. 
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which is a sub-bundle of E(M) . Locally, ECu, U) looks like Uy Xe, 


1.5 Vector Fields, Connections, and Lie Derivatives. 


In this section we give a fairly detailed discussion of vector 
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tangent bundle of the bundle of linear frames. Having defined vectors 
everywhere on G and T(E(B)) by means of left-invariant and fundamental 
fields, we talk about connections between fibres of E(B) , and use the 
concept of the Lie Derivative to introduce covariant differentiation of 


tensors on T(M) 


In dealing with the bundle of frames, it is helpful to keep the 


diagram below in mind. 
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For generality, we take G=GL(3) and assume that G is 


coordinatized, i.e. aeG _ is represented by CAND) fo ues 
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Now if X| is a tangent vector at the identity of G , we may construct a 
vector field X which is left-invariant, by defining the generic element 
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Conversely, if X is left-invariant, then it must be of this form. The set 
of left-invariant vector fields has as many elements as the dimensions of 
G_. In particular, the set of left-invariant vector fields on GL(n) 
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The set of left-invariant vector fields on G_ forms a Lie- 


Algebra, g , where the bracket operation is given by: 
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where x is a vector field on the fibre E(p) . Since coordinate changes 
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1.6 Material and Riemannian Connections. 


The discussion above on connections may be easily specialized to 
define material connections. We restrict our bundle of frames E(B) to 
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is a Lie subalgebra of g. This induces a Lie subalgebra ig(U) mon EUb.U) a, 
with corresponding fundamental fields which are restrictions of the fields 


considered above. 


Suppose we have a curve X(t) through p on M. Let (uv) 
be a coordinate chart around p_ and let the coordinates under Yo be 
{x'} . If Xe) is a smooth natural lift of A(t) , then the tangent vector 


TOME ye ah u = (p,u) (ns E(B,U). sis, cf. (125.50). 
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Now this vector is in H« for some connection [I iff 
pil =P ieee ae (1.6.2) 


We then say the frame U is moved along A(t) by parallel transport. 
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Parallel transport of vectors induces isomorphisms of the tangent 
spaces of M along A(t) . If the induced isomorphisms are material 


isomorphisms, we say [ jis a material connection. 


Now suppose we have a material atlas on M , and a smooth cross- 


section oO on Cu by » given by 
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isotropy group, we see that [ is a material connection iff 
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Associated with a given connection are two further tensors, the 
curvature tensor 
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and the torsion tensor 
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A given material connection is (locally) flat if for every pe M 
there is a local coordinate chart in which the connection symbols vanish. 
It is known that a necessary and sufficient condition for a connection to be 


flat is that 
Faull) ay ae 8) Cis On) 


There is no reason why a given material connection should have either property. 


Referring to the definition of local homogeneity, we see that if a 
body has a flat material connection, it is locally homogeneous, since there 
exists a coordinate chart around each point, in which the material and natural 
frames coincide. For solid bodies, we also have the converse, i.e., local 


homogeneity implies flatness. 
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Suppose we have a solid body with a reference atlas U_ relative 
to which the isotropy group is a subgroup of S0(3) . U_ is called an 


undistorted atlas. We may introduce a Riemannian metric on U _ by taking 
G CH, Vor Cin aes ie (v) Cleon iy 
p O,P O,P 


where. u,v € ue : = ry 6P? is a reference chart, and the Inner product 


. 

a,P 
is the ordinary Euclidean inner product on R .. The orthogonality of the 
isotropy group shows that this is an intrinsic metric at p . Then we can 
extend Cu, to an intrinsic field on U . Suppose we fix this atlas. 


Given any material connection, it follows that the covariant derivative of 


S| y with respect to the material connection must vanish. We get 
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where the I's are material connection symbols. If the material connection 


is torsion-free, i.e., symmetric, then it is uniquely determined by G. 


On the other hand, if we introduce an arbitrary intrinsic Rieman- 
nian metric g , then we define the Riemannian connection to be the unique 


torsion-free connection such that 
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where in this case the I's are the ordinary Christoffel symbols. If This 
Riemannian connection is not material, then no material connection is 
torsion-free. If one of the material connections is torsion free, it must 


coincide with the Riemannian connection. 
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Looking at the example we gave earlier of a transversely isotropic 
spherical shell, we take a configuration such that the intrinsic Riemannian 
metric coincides with the Euclidean metric. The induced connection is not 
material, so there is no torsion-free material connection. Since any locally 
homogeneous body has a torsion-free material connection, this is an example 


of a body which is not locally homogeneous. 


There may be several undistorted reference atlases for a given 
solid body, and so more than one intrinsic Riemannian metric can be induced 
in the above manner. Coleman and Noll have shown that if the solid is 
isotropic, the metric is unique to within a multiplicative constant. Since 
this constant will not affect the Christoffel symbols, a solid isotropic 
body has a unique Riemannian connection. Since this connection maps ortho-~ 


gonal frames to orthogonal frames, it is a material connection. 


Generally the connection will not be flat. If the body is also 
locally homogeneous, then it is. Given any locally homogeneous solid 


body, isotropic or not, there is an undistorted atlas U = f }. By 
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local homogeneity the field of reference maps ry corresponds to the 
induced local configurations of the coordinate map vy One Us The compo- 
nents of the intrinsic metric relative to by are On » so the connection 


is flat. As pointed out above it is also material. 


1.7 Laminated Bodies. 
An isotropic solid body B is called a laminated body if it is 
a disjoint union of a collection of two-dimensional submanifolds Le > cal led 


the laminae, 
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where Bis an index set, and where each Le has the property that in 
some neighbourhood of it, there exist configurations by whose induced local 
configurations carry the intrinsic metric G on L, onto the Euclidean 
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metric, i.e., 
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different set of v's may have to be chosen for each neighbouring lamina. 


The v's are called initial configurations for Le . 


Physically, we may think of the body as made up of a set of 
infinitely thin, "pre-stressed" sheets, or of a set of sheets deformed 


separately and then put together. 


The body is inhomogeneous. There is in general no global config- 
uration (deformation) of the body which will return it to an undeformed 


state. 


It is shown in [7] that for certain classes of laminated bodies, 
+here aretniversal solutions (possible deformations), which are the same as 
those for homogeneous bodies. One such class, and a universal solution for 


it, is given below. 


A laminated cylinder is made up of thin cylindrical shells. We 
assume that there exists some set of initial configurations such that rela- 


tive to a cylindrical coordinate system (R,0,Z) the components of the 
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intrinsic metric G (or g, to bring the notation into line with chapter 


2), form the matrix 


g 0 0 
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5 a2 =33 


where the g's are functions only of R. We deal with incompressible 


bodies, so 


det (g"")) = 1/R* Grages) 


From (1.3.4) we have a representation of thestress for an isotropic 


elastic solid. We may write it as 


T= ~=pl + ¢,8 + f_,B Cird* 2) 


where I is the metric tensor of the space and B_ is the deformation tensor 
relative to some fixed undistorted reference configuration. If we choose 
this configuration to be the one in which the intrinsic metric and the spa- 


tial metric coincide, then we have 


t= -pl+fg+ f(g). (1.7.6) 


We may write this simply as 


t=-pl+S , Ghee 
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and the S's are functions only of R. 
Suppose we deform the body according to 


mee dialer, Lis yestivar lye ois ge he (1.7.9) 
where A,B,C,D,E,F are constants satisfying 
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In the deformed coordinate system (r,8,z) the intrinsic metric q has the 
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so we have incompressibility satisfied. 


The equilibrium equations for a deformation are 
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Since we are dealing with an isotropic solid, the Riemannian connection is 
the material connection we need for covariant differentiation. Several of 
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We see that this deformation is in equilibrium under a pressure given by 
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CHAPTER 2 


The Equations of Equilibrium for a Perturbation 


Acc Deformations. 


We deal exclusively with simple, materially uniform, solid, 
isotropic bodies. Suppose we have an initial configuration X fora 


body B , and a deformed configuration x where 


vey ae (2.1.1) 


both coordinate systems being rectangular Cartesian. The deformation 


gradient F has representation 


ax! 
F= (25), (2.162) 
xX 


and the left Cauchy-Green tensor B has components 
Ss 
BY = — 6 C2 oe 


referred to the coordinates {x'} 


In order to determine a standard way of examining small deforma- 
tions on large, we bring this notation into line with that of Green and 
Zerna [2]. We must then consider general curvilinear coordinates. Let 


{o"} denote such a set in the initial system with base vectors 
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Under the deformation map from X to x the induced map of vectors gives 


a set of base vectors 


a 
g = 8 (2.1.6) 
06° 9 
and a corresponding metric 
Aes 
es oe om (2: 17) 
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A state of deformation may already exist in the initial configura- 
tion, and some care is required in considering the composition of deformations. 
lf we have three configurations Xo » X , and x and deformation gradients 


as shown diagrammatically below, Then 


Ff > F °F : RG bes 8. 
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Let B. denote the Cauchy-Green tensor in x due to the 
deformation from Xo to x , and By the tensor in X due to the 


deformation from Xo to X. Then 


veg 
pre = 2X Be (2.1.9a) 
dX OX 
A B 
mit _ dX» 0X -| 
(B. Ra = is NEE ) AB C2 Ue a) 
ax, Ox 
Ge Mesh. Ss 
(Bays = Oe ee eg (2.1.9¢) 
9X 3X 8X ca . 
To follow Green and Zerna, By will be given in terms of a 
coordinate system {6'} and B in terms of a system {61}. In the case 


of a laminated solid the tensor By coincides with the induced metric g 
and we assume this is known in the initial state. To be definite we take 

the components of g_ to be gt relative to a given curvilinear coordinate 
system cr in the intial state. Then relative to the {6°} in the 


deformed state, g has components ae given by 


eee Eric: 
ge. Beg (2.1.10) 
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Let Bo have components Be relative to the convected coordi- 


nate system {0%} in the deformed state. If we set {e"} = {6} then 


some computation yields 


Bore gue (2.1.1 1a) 
@') = gw) (2.1. 11b) 
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(Bo = 6 aoe Cae (2.1.1 10) 
p 
and we find the invariants 
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I, = det B = G/g : (2.1.12c) 


where G,g are the determinants of ((6%)) ; Ca a3 respectively. 


2.2 Small Deformations on Large. 


The deformed body is set in a curvilinear coordinate system {oe} . 
Suppose we impose an additional small deformation. We write the new position 


vector of a point p as 
Borage wits se, C2. by 


where ris the former position vector, e¢ w(r) is the deformation field, 
and e« is asmall parameter. We derive the first order theory. In general, 


we denote a perturbed quantity Q. by 


T=Q+eQ (a8) 
The new tangent vectors are 
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We may now determine the equilibrium equations. Let as be the 


components of the stress tensor in the {0%} system, and let 


qo = 798 eG Gp ; (2.2.12) 


Then 
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The A's are not symmetrical in general. 
The Cauchy condition for the perturbed state is given by 
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Since the body is an isotropic solid, the material and Riemannian connec- 


tions coincide, so we may use the symmetry property of the connection, and 


the fact that 
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and the equilibrium equations become 
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The boundary conditions also have a simple form 
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We now digress, to study the form of the t's and A's in 
greater detail. For an isotropic solid we have, using (1.3.4) and (2.1.11), 
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In particular, (see Appendix |) we have a strain energy function W = and 
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We may determine the perturbed form of this equation fairly easily. Consid- 
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The constitutive equation for a Mooney-Rivlin material is of the 
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ILS A Particular Deformation. 


We consider a laminated right-circular cylinder under vertical 
compression. If we choose cylindrical coordinates (r,8,z) in the deformed 
configuration, and let {06°} be (R,0,Z) , we have a universal solution 


of the form (1.7.9), namely 
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The deformation metric gor then takes the form (cf. (1.7.11)) 
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As shown in section |.7, the incompressibility condition is satisfied, i.e. 
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As shown in section |.7, this deformation is in equilibrium under a pressure 
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We write the integrand as 
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We wish to impose a perturbation on this compression, but the 
general system (2.3.10) gives perturbation equations of unreasonable complex- 


ity. We make some simplifying assumptions first. We assume 
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We also diagonalize the intrinsic metric by setting 
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WIth this choice of the metric, we see that the T's are diagonalized. If 


we put (2.3.13) into (2.3.9), we get 


COi=.0'°, (Zea) 
whatever The values of © and Y¥ , and hence ~ is the only non-vanish- 
ing component 

= oa agi) (2.3.15) 
Cg) -s) 
We let 
Gok gt | | GS ie 


BBeESUNS Ltt eubimte sooo 2.5. wr and Cl.7.lD) TO obtain 


Cin = Mero i) PASM} 
iF 
Daze Out. C208) 
In the case of a homogeneous body, 
heel ie ae (2.3.19) 


Suppose now that a small perturbing force due to gravity acts on 


the cylinder. It has components 
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We assume that the perturbation has vertical and radial effects, but no 


angular effect, so that the deformation field is given by the vector 
we) = (u,0,v), (2.3.21) 


where u and v_ are functions of r and z. We can now apply the pertur- 


bation theory of the last section to obtain the new equilibrium equations. 
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We also have the incompressibility condition (2.2.37): 
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lf we take the partial derivatives of (2.3.29) with respect to r and z, 
the resulting identities may be substituted into (2.3.27) and (2.3.28). The 


set of equations then becomes: 
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These equations, together with appropriate boundary conditions, 
determine the perturbations u,v uniquely in terms of f , although not in 


a closed form. 
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CHAPTER 3 


Compression of a Laminated Cylinder Under Gravity 


Bel Description. 


Consider a vertical right-circular cylinder of initial height 
hy » radius ee supported at its horizontal faces by planes. We assume 
that no friction exists at these faces. The cylinder is laminated in such 


a way that the intrinsic metric coincides with the one GiVens I Ne teeoeho): 


We first suppose the cylinder to be compressed to height ab 
with the ends maintained in a horizontal position. Essentially this problem 
was solved in the previous chapter. The stress system is given by (2.3.10). 
Our purpose here is to investigate the effect of a small perturbation on this 


solution. 


We refer to a paper of Vaughan [5] in which a similar investiga- 
tion was carried out for a homogeneous cylinder under gravity. The basic 
assumption is that gravity effects are small and can be treated by first- 
order perturbation theory. We may compare our results with those of Vaughan 
[5] by setting our intrinsic metric equal to the ordinary spatial metric. 

As we shall see, our solutions reduce to those in [5] if we deal with the 


homogeneous case. 


We review some of the results derived in chapter 2. The compres- 


sion is given by (cf. (2.3. 13)): 
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The intrinsic metric is 
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and we define f , a positive function of r , by 
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lf g is the constant perturbing gravity force, and u,v are 
the horizontal and vertical perturbations, respectively, then the perturba- 


tion is in equilibrium under a pressure field 5 given by 


dp _ , df du, > 9 du u , % du 
spe NP apes pes eee ees NUH pe ite he ces (Z.5.0a) 
G fmnoz 
a 2 : 
ejay a ahs, Che See a) OV dv _ 
Sn Vege oo Cah aye ar epee 2: Pg ° (275.50b) 
The incompressibility condition is given by 
OUR, Ue, ON GE Ey : (225-500) 
eo lds Zz 


The boundary conditions are given by (2.2.30). There are no forces acting 


on the curved face r = ro» SO using (2.3.26) we have 
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The ends of the cylinder are frictionless and are maintained horizontal by 


the action of the planes, so we have 
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The problem then is to solve (2.3.30) for pB,u , and v , subject 
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where the integration is carried out with respect to z . The equilibrium 


equations (2.3.30a,b) now become 
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We now derive a compatibility condition which must be satisfied 
if the equations (3.2.5) are to be self-consistent. Integrating (3.2.5b) 
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The compatibility condition is obtained from (3.2.7) and (3.2.5a). After 
some reduction, we find that we must have 
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We now study the possible forms that the function Z may assume, 


and reduce this equation to an ordinary differential equation for T . 
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If Z is a polynomial of degree v , then | | Z is of degree 
v+2 . To satisfy (3.2.8), the two highest powers of z must vanish, there- 
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The only solutions of (3.3.1) which satisfy (3.1.4) are 
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where a and 8 are arbitrary constants. 
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using (5.5.2), and since n is arbitrary, this condition can be met. 


Therefore we have as possible deformations 
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where a,8 are arbitrary constants. In the homogeneous case where f is 


a constant, these deformations reduce to 
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which agree with those of Vaughan L5]. (To satisfy (3.1.5), set B=ye=0O). 
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The equation (3.4.4) may be written 
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in the homogeneous case, and since the operators commute, we obtain two 


linearly independent non-logarithmic solutions from (3.4.7), namely 
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where Q, and a, are constants. This form leads to Vaughan's solution. 
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boundary conditions may be satisfied. 
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3.6 The Power-Series Approach. 


We were unable to determine the solutions of (3.4.4) in closed 
form. Atempts to derive explicit expressions for T by means of power- 
series, by expanding f in an arbitrary power-series, also achieved no 
success. We thus decided to investigate solutions corresponding to an f 
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qT, » corresponding to (3.4.9b), is obtained by choosing 
area, a, = E ; (356.20) 


At the origin, the inhomogeneous and corresponding homogeneous solutions match 
in value, slope, and concavity. The boundary condition (3.1.4) is thus 


satisfied, as indicated in section 3.5, 


A program was written to compute the first fifty ays for various 
values of L and M. There was some oscillation in signs of the a's for 


larger values of M , but their absolute values decreased extremely rapidly, 
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and a further division by m! to obtain the coefficients of s” in (3.6.10) 


ensured that the last thirty or more a's had negligible effect on solutions. 


The power-series method was an inefficient use of computer time. 
It was only used for two specific purposes. I+ was used to calculate solu- 
tions and their first three derivatives at a small s-value. These results 
were then used an input for a faster and more efficient program. The method 


was also used to calculate the two solutions corresponding to a particular 


M- value, M= - a opt which caused the function oO) TO Vanish shathis 


particular result was used in studying the behaviour of the determinant of 


The system (3.6.9). 


3./ The Numerical Solution. 


The equation (3.6.4) was converted to a system of four first-order 


ordinary differential equations. Setting 


Casts Mest, C2 oh asl 


and dividing (3.6.4) by s*o » we are led to the system 


4 
* 
) 
>) 
° 
4 
Pen 


ie 0 0 | 0 aR 
Vite 0 0 0 | Te 

| L> M 24 | L 2 
KKK = aa mes os pe a i ees KKK 
h Ze BES eA ae! a ma Eehns me 2M] - —[l+2Ms] t 


JER Eys | 


(Ol,8.0) al “2 to etaetolt team Gh Oi or Ga ye. OMEN TR Ge oe brs, : 


ENO THIOS NO TosttS Gti gi igan Gan eS oa yt feet ont tent beuens ; 


» 


La 
mu rat ; saul Tr On; WO @iw bo ‘Ter see eweq- AT ; 
; ? 
Pn von 
“eae Se of bequ 26w (| /eiMwn ol*iseqe awh 36) boey ving Sie ine 
‘ 7 ( 
ws : ‘Mi oves ' 6 TD tevitévined aay? fart? +hed? bos ange 


iSivtS eranr Gris. +o wor Tugo! ea beev Aent a9W - 


7 


aa ne on Get Gn! etéluolts of beau cele. Sam 


b = 
3 Holi, gee ~ eM ey see 
u ) Lal 
2 ¢ ; . ie 
ivnniee 92%) ontybu' bed. Sew tivee? “elu Tega 
~ < 
(@.3,8) mataye = - 
all 
tN int 
Moltule? ltolvemyetT, Fae 
Pe LSS Se : 
7 6 
, ' 08 fod 2aw'¢ iE) me I supe ewiT ' 


nh pee lohtpups bei teaneltio yientbre 


yo (PB, 3.8) enlbivib bnay 


‘- 


7 ee 


This system was solved on an IBM 360-67 installation using the 
"DASCRU" program from the IMSL library. The program uses a modified Runge- 
Kutta algorithm. It is the fastest program available for this purpose and 
has an accuracy of approximately machine-word-length. Required input 
includes the order of the system, an initial-value vector, and the function 
which determines the final row of the matrix. Output is a solution vector 


at a specified s - value. 


The "DASCRU" program could not handle the system (3.7.2) when s 
or 9 vanished. This is why the initial-value vector was calculated by 


power-series. 


In order to evaluate boundary conditions, it was necessary to choose 


a specific value of Ce We chose 
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In terms of s , the interval over which solutions had to be calculated 


became, for each nn, Wezauhes - The constant S| in (3.6.7) became 2.5. 


For each choice of M and L _, the power-series method was used 
to calculate an intial - value vector for each linearly independent solution, 
for each value of n . Each initial-value vector was used as input to 
"DASCRU", which gave solutions at ten equal intervals outward from .25n7 


To 2.5n° , using the output from each step as input for the next. The 


solutions are tabulated in Appendix 2. 


Appendix 2 also includes various intermediate results, mentioned 


below, leading to the final form of the perturbations. 
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The solutions and their first two derivatives were positive 
everywhere, in fact the second derivative had positive minimum. Thus there 
is a strong indication that the general equation (3.4.4) is non-oscillatory. 
The solutions for negative M are apparently bounded below by the homogen- 
eous solutions (the Bessel functions), while the solutions for positive M 


are bounded above by them. 


The boundary values were used as input for a program designed to 
solve the system (3.6.7), as indicated in section 3.5. For all M - values 
for which 0 did not vanish, the determinant was positive and became very 
large for the higher harmonics. We were thus able to obtain unique solutions 
of. the form (3.5.11) for the significant harmonics. The solutions for 
Ce O were calculated by power-series. The determinant of the system was 
then found to be a three-digit number representing the difference of two ten- 
digit numbers. Therefore, within the limits of machine accuracy, the deter- 
minant of the system (3.6.9) vanished. We were thus unable to obtain a 


Ebimuring!s solution of athe form (5.5.11). 


The unique solutions obtained above were then converted into a 
form suitable for substitution in (3.4.1). The series contributions to the 
perturbations were then calculated, to three or four harmonics, i.e. to 
within 1.5% . The final results included various constant scaling factors, 


but comparisons could be made. 


The horizontal perturbations were up to 45% smaller for M = -2 , 
than for the homogeneous case, while those corresponding to M= +.2 were 
up to 45% greater. There were also significant, though not monotonic, 


effects on the vertical components. In general, a negative M - value gave 
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less vertical perturbation as well. 


The polynomial terms in (3.4.1) represent the first-order form 
of a universal solution of the type (2.3.1). Hence it may be assumed that 
these terms are included in the initial compression, so that the actual 


perturbation is obtained by setting a= 0 in (3.4.1). 


We have thus shown that a laminated metric of the type (2.3.13) 
prohibits any other deformation of polynomial type from being superimposed 


on a solution of the type (2.3.1), that is, a simple compression. 


We have also shown that the laminations have a considerable 


effect on the magnitude of the perturbations due to gravity. 
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APPENDIX | 


Thermodynamics and Symmetry of the Constitutive Equation 


This is an outline of the work of Rivlin [3] and Ericksen and 
Rivlin Ll], in which they developed the simple form of the constitutive 


equation (1.3.4) for an isotropic solid under an isothermal deformation. 


The material time-derivative of the position vector x is denoted 


by v = x! We deal with rectangular Cartesian systems. 


Denote absolute temperature by K , surface forces by t' , body 
forces by f' » the rate of heat generation by r , heat flux by q, 
original density by Po» internal energy by U , and entropy by S_, all 


quantities defined per unit mass or per unit undeformed area. 


The balance-of-energy equation then becomes 


(tiv, -q)dA = | p (UHV v dV ’ (Al. 1) 


| p (rtf v.)dV A | 
V V 


A 
while the Clausius-Duhem inequality is given by 


; 0, 
| p, 8 av > | =r dv - | zaA (Al .2) 
V V A 


lf in an elementary tetrahedron we let nih be the force vector 


acting on the surface whose normal is My > and a be the heat flux through 
that surface, then in the limiting case as the volume decreases to zero, we 


must have 
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lf we impose a constant translational velocity c , we get 


i i peel Lk iA eA 
tv, ot Cha Gar ny. + 1 NAc; qn, (Al .4) 
and therefore we must have 
aS wy (Al .5a) 
PSA 
See Se RAL.) : (Al .5b) 


Putting these results into (Al.1), we obtain, after some reduction, 


A 
[A 


| iA A . 
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Again considering the effect of an imposed velocity c , we conclude 


iA = 


pf +m = ply Cale Ta) 


iA AVS Seca ae 
T Vila ROE CC aed [A ~ POU . (Al. 7b) 


Using the divergence theorem, (Al.2) becomes 


A A 
p S ? ROME : lA + é SA (Al.8) 
o — K K Ke : . 
and using (Al.7b), we find 
e e iA aK, 
PCU-KS) <7 Vila - r (AI.9) 


Define the Helmholtz free-energy by 


A =U - KS's (Al. 10) 


= B28% 


Since the deformation is isothermal, K vanishes, and in (AI.9) we have 


(Al. 11) 


The energy A_ depends explicity only on the deformation gradients, and not 


on the time, so (Al.11) becomes 


A 

OW Le We! J eas 

°F T 6, y iS ac Ky, (Al. 12) 
a 5p 


To satisfy this identically for all deformations, we must have 


aia ees (Al. 13a) 
eee ee ee (Al. 13b) 
© J ij 
"A 
The stress tensor t'! is related to the t's by 
iA Jj 
Sp Oe Mae ae 35 
rid = (AL. 14) 
det F 
where F is the deformation gradient matrix. Thus 
Juba eal ec (AI. 15) 


»A ax”, 
where 09 is density in the deformed state. 


By the principle of material frame indifference, the energy A 
must be invariant under all orthogonal changes of frame. It may then be 


shown that A must depend explicitly only the the C defined by 


PO. 


G = x! 


PO (AL. 16) 
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Equation (Al.15) then reduces to 


do 8A 
Q CE 


Bie i J i 
x (x 9 x" 55 To »p x (Al.17) 
This is a standard result for all bodies. It may be simplified 


further by noting that tT is invariant under transformations of the isotropy 


group. 


lf the body is transversely isotropic, an element of the group is 


| 0 0 
((HE)) =! 0 cosw  sinw ; (Al. 18) 
0 -sin w COS W 
and if we define 
~ : Q 
Cor Cog He (Al. 19) 


then we must have 


pp = Sawer aa = : (Al .20) 


This identity can only hold tf there are certain relationships 
among the components of C , so that A_ depends explicitly only on certain 
functions of the components of C . These functions also satisfy certain 
identities involving powers of C and powers of B , where 
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In particular, for an isotropic material, the group element in 
(Al.18) is any member of the group S0O(3) , and the identity (AI.20) 
requires that A depend only on the three invariants of B , so that the 


final form of (Al.1I7) becomes 


TE oeCHaY +i Hebe! op!, (Al.22) 


where the n's_ are functions of the three invariants. Since Bis non-singu- 


lar and satisfies its characteristic polynomial, 
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another form of the constitutive equation is 


T= EL+EBHE, Be (Al .24) 
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APPENDIX 2 


Numerical Tables 


In this appendix we give a sample of the numerical results which 
illustrate the methods and conclusions of Chapter 3. In particular, we 
give horizontal and vertical components of the displacements of representa- 
tive points in a laminated, compressed cylinder perturbed by gravity. The 
quantities calculated varied smoothly with respect to all parameters, and a 


condensed sample is sufficient to indicate the behaviour of solutions. 


The results given are the series contributions To the displace- 
ments (3.4.1), up to the fourth harmonic. The results are also scaled, and 


are thus related to the u,v of (3.4.1) by 


teh 

ee = (u-ar) 
given 80g 
on h, 

= (v+2az) 


Ygiven  ~8pg 


The representative points lie on two concentric cylindrical shells, 
v10 


v3 ; 
an interior one of radius a and the outer surface, radius srorgcts i For 


each harmonic, the corresponding s - values are 75% and 2. on 


This appendix is in three sections, each one corresponding to a 
particular lamination characterized by setting M S52) 5 Up eut coe 
(3.6.6). The second section, M = 0 , represents the homogeneous case. The 
constant L_ in (3.6.5) is assumed to be ./29 , representing a compression 
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In each section, the first page of tables gives the solutions qT 


and qT, of (3.7.2) for each harmonic. For the case M = 0 , the required 


solutions and derivatives were obtained from Bessel-function subroutines, 


which were simpler than DASCRU. 


The second page includes the determinant A of the system (3.6.7), 


and the constants determined by it, for each harmonic. The unique solution 


(3.5.11) thus obtained is multiplied by né and listed, along with its 


derivative. These are the actual coefficients of the sine and cosine terms, 


80g 


on'h, 


except for the scaling factor 


Finally, the horizontal and vertical components of displacement 
are given, for various z - values. A diagram showing these displacements 


is also included. 
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Output of Bessel Function Subroutines 
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